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COMPARISON OF CONVENTIONAL AND MICROASPERITY ELASTOHYDRODYNAMIC 


LUBRICATION OF A BALL SPINNING IN A NONCONFORMING GROOVE 
by Charles W. Allen,* Dennis P. Townsend, and Erwin V. Zaretsky 

Lewis Research Center 

SUMMARY 

The microasperity elastohydrodynamic lubrication analysis for a ball spinning on a 
flat plate was extended to the elliptical contact case of a ball spinning in a nonconforming 
groove. The microasperity analysis was compared to the conventional elastohydrody- 
namic analysis of a spinning ball in a nonconforming groove. Rheological models for a 
di-2-ethylhexyl sebacate, a super -refined naphthenic mineral oil, and a polyphenyl 
ether (5P4E) were constructed, by using both analyses, from spinning torque data ob- 
tained in a previous investigation. The value of the lubricant pressure-viscosity coeffi- 
cient that makes the data fit the analyses of the fluids differs somewhat from published 
values. For all three lubricants, an exponential composite model best described the 
lubricant rheology. Good agreement existed with the experimental and analytical values 
of torque for both the conventional and microasperity elastohydrodynamic analysis for a 
spinning ball in a nonconforming groove. The polyphenyl ether lubricant exhibited a 
negative secondary slope in the composite exponential model, which would imply a de- 
crease in viscosity with pressure above 117x10 newtons per square meter (17 000 psi). 


INTRODUCTION 

In recent years, a complete evolution has occurred in ball bearing design and anal- 
ysis. Early work assumed that at the ball-race contact, spinning and rolling would 
occur at one raceway and only pure rolling at the other raceway (refs. 1 and 2). This 
work was based on the premise that there was a constant coefficient of friction at the 
ball-race contact. Work reported in reference 3 and subsequently in references 4 and 5 
indicated that the friction was not dependent upon a single coefficient of friction but var- 
ied according to the contact stress, contact geometry, and lubricant type and viscosity. 
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As a result, no simple coefficient of friction could be applied to a ball-race contact. 
Reference 6 recognizes this phenomenon and presents an analysis for the kinematics of 
a ball bearing considering elastohydrodynamic (EHD) effects and assumed lubricant 
rheological properties. However, these properties are not completely defined, resulting 
in some second-order inaccuracies in the analysis. Reference 7 develops a theoretical 
analysis based upon the elastohydrodynamic theory of lubrication for spinning torque. 
This model was verified experimentally. Subsequently, a microasperity analysis for 
elastohydrodynamic lubrication of a ball spinning on a flat surface was developed 
(ref. 8). By utilizing these two analyses, it is possible, based upon experimental mea- 
surements, to derive a rheological model for a test lubricant. 

The work reported herein was conducted (1) to extend the microasperity analysis 
for elastohydrodynamic lubrication of a ball spinning in a nonconforming groove; and 
(2) to determine the pressure -viscosity models for a di-2-ethylhexyl sebacate, a super- 
refined naphthenic mineral oil, and a polyphenyl ether by using a conventional EHD anal- 
ysis and a microasperity EHD analysis. The spinning torque data used were obtained in 
the NASA spinning torque apparatus and originally reported in reference 3. These data 
were verified in the modified spinning torque apparatus described initially in refer- 
ence 4. 


APPARATUS AND SPECIMENS 
Spinning Torque Apparatus 

A spinning torque apparatus (see figs. 1(a) and (b)) as reported in references 3 and 4 
was used for the tests reported herein. The apparatus essentially consists of a turbine 
drive, a pneumatic load device, an upper and lower test specimen, a lower test-housing 
assembly incorporating a hydrostatic air -bearing, and a torque-measuring system. In 
operation, the upper test specimen is pneumatically loaded against the lower test speci- 
men through the drive shaft. As the drive shaft is rotated, the upper test specimen 
spins in the groove of the lower test specimen. This causes an angular deflection of the 
lower test-specimen housing. This angular movement is sensed optically by the torque- 
measuring system and is converted into a torque value. During a test, the torque is 
continuously recorded on a strip chart. 
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Turbine drive air 


/- Drive turbine 



(a) General cutaway view. 

Figure l - Spinning-torque apparatus. 
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Front view Isometric view 
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(c) Lower test specimen. 
Figure 1. - Concluded. 


Specimens 

The upper test specimen is a conventional 12. 7-millimeter- (1/2-in. -) diameter 

bearing ball made of SAE 52100 steel having a nominal Rockwell C hardness of 61 and a 

-8 

surface finish of 5X10 meters (2 juin.) rms. The lower test specimen (fig. 1(c)) is a 
12. 7-millimeter- (1/2-in. -) diameter ball from the same heat of material as the upper 
test specimen, which is modified by grinding a flat on one side and a cylindrical groove 
of radius Rq (fig. 1(c)) on the other. The groove simulates the race groove of a bearing. 
The axis of the groove is parallel to the flat. The groove radius expressed as a per- 
centage of the upper-ball diameter is defined as the ball-race conformity. The speci- 
mens used in these tests were ground to ball-race conformities of 51 percent. The 

-8 8 

surface finish of the cylindrical groove was approximately 5X10 to 15X10 meters 
(2 to 6 ju in. ) rms. 


THEORETICAL ANALYSIS 

An elastohydrodynamic analysis for the determination of the torque of a ball spin- 
ning in a nonconforming groove, shown in figure 1, was developed in reference 7. In 
this analysis, the area of contact is elliptical, as shown in figure 2(a). In this ellipse an 
elastohydrodynamic film can theoretically be predicted, due to the spinning velocity, ex- 
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(a) Contact ellipse for ball in nonconforming groove. 


— *-j |— dy 



(b) Elemental roller. 


Figure 2. - Ground lower ball showing contact ellipse and elemental roller. 
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cept within the inscribed circle of radius b. The frictional force over the ellipse is de- 
termined by making the system into a number of elemental rollers, as shown in fig- 
ure 2(b), and determining the frictional force on each roller. By using this method the 
frictional force dF in an elemental roller of width dy sliding at a velocity of yo> is 

dF = 2b^y dy (!) 

h 


where the equivalent viscosity p as given in reference 9 is 




b’ 



p dx 


( 2 ) 


The film thickness h from reference 10 is 


h = 


1.6o°- 6 E’°- 03 R 0 - 43 ^o Wy 


,0. 7 


W 


0.13 


(3) 


where 


E’ - 


a 


For steel, E’ is 22.3X10 10 N/m 2 (32.3X10 6 psi). 

The assumption of a Hertzian stress distribution leads to a load distribution W. 


W - 



(4) 


where P is the normal load. 

The lubricant viscosity on the surface at location x is dependent upon the contact 
pressure at that point. If a Hertzian distribution is assumed, the contact pressure is 
given by 


s = MP 
7rab 



(5) 


6 



The actual viscosity depends upon the pres sure -viscosity function. 

The friction moment M.^ about the Z-axis due to the section of the contact ellipse 
for y > +b and y < -b is obtained by integrating the product of the friction of the ele- 
mental roller and the moment arm; thus 


M 


1 



( 6 ) 


A numerical solution is used for this equation since a closed -form solution does not exist. 

There remains the area within the inscribed circle of radius b where a true elasto- 
hydrodynamic film theoretically is impossible to maintain. Therefore, the film thick- 
ness h is assumed to be the same as that at y = b. If the inscribed circle is then divi- 
ded into elemental rings of radius r and width dr, the torque Mg for the inscribed 
circle becomes 



where is the equivalent viscosity at y = b. 

M g = M 1 + Mg (8) 


The preceding paragraphs outline the general approach used for the conventional 
elastohydrodynamic analysis. In using this analysis for the friction torque of a ball 
spinning in a nonconforming groove, it was found in reference 7 that the commonly used 
exponential law p = fJ- 0 e01 ^ for the pressure-viscosity relation gave torque values that 
were much higher than experimental results. However, a good correlation between 
theory and experiment was achieved by using a composite pressure-viscosity relation 


,o;S 




^ 


(u. = fz Q e u where S < 
o?Sj+/3(S-S^) 


where S > S. 


(9) 


Equation (8) in the preceding analysis cannot be used for a ball spinning on a flat 
plate because of the absence of a converging film; that is, two flat surfaces theoretically 
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cannot form a hydrodynamic film. However, the lubrication of a ball on a flat plate has 
been explained by a microasperity EHD theory (ref. 8). In this theory an EHD film is 
formed between the asperities on the two flat surfaces, as shown in figure 3. For an 
individual asperity, as shown in figure 4, the film thickness is given by 

h = 0.67 (ap Q v) 2//3 p 1/3 (10) 

Considering all the asperities contained within the contact ellipse the total area of con- 
tact is given by 



where is the average Hertz pressure over a typical asperity. If the fluid is as- 
sumed to behave in a Newtonian manner and the velocity gradient is assumed to be linear, 
the shear stress r in the fluid at the top of the asperity is given by 



(a) Two rough surfaces in contact. 



(b) Equivalent rough surface on plane. 


Figure 3. - Surfaces in contact. 


8 




Figure 4. - Mode of asperity lubrication. 


Assuming (1) that the fluid between the asperities does not contribute a significant 
torque and (2) that the fluid viscosity is constant over each asperity peak, the total spin- 
ning torque imposed by the microasperity film-over the whole contact ellipse is given by 



This may be integrated by using numerical methods. 

From the preceding discussion it can be seen that the spinning torque for a ball spin- 
ning in a nonconforming groove may be determined by either of two methods. The first 
is conventional EHD analysis using equation (8), which assumes a constant film thickness 
over the inscribed circle. The second is microasperity EHD analysis using the method 
developed for a ball spinning on a flat plate in reference 8 and extended to a ball spinning 
in a nonconforming groove. 

In addition to the effect of the lubricant within the Hertzian contact region, the effect 
of the lubricant outside of this region must also be considered. The expression for this 
moment is given in reference 5 as 



This may be integrated numerically over the region outside the contact ellipse. 
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In reality, the conventional EHD and microasperity effects are superimposed but, 
since the coupling between the two systems is extremely complex, the conventional EHD 
and microasperity EHD effects will be examined individually. The effect of the lubricant 
lying outside the contact zone will be added to each. 


RESULTS AND DISCUSSION 

The spinning torque reported in reference 3 for a di-2-ethylhexyl sebacate, a super - 
refined naphthenic mineral oil, and a polyphenyl ether was examined by both the conven- 
tional and microasperity theories for a ball spinning in a nonconforming groove. For the 
di-2-ethylhexyl sebacate, good agreement was obtained between the experimental results 
and each theory by use of the following parameters: 

a - 2.03X10 -8 m 2 /N (1.4X10 -4 (psi) -1 ) 

j3 = 1.45X10' 10 m 2 /N (l.OxlO -6 (psi) -1 ) 

Sj = 275x10® N/m 2 (40 000 psi) 

fi = 16X10 -3 N-sec/m 2 (2.3X10 -6 lb-sec/in. 2 ) 

The pressure -viscosity relation given by the preceding values is portrayed in figure 5. 
The spinning torques for the conventional and microasperity EHD analysis are shown in 
figure 6. The value of a used here is higher than the published values of 

a = 0.58X10 -8 m 2 /N (0.4xl0 -4 (psi) -1 ) (ref. 11) 


and 

a = 1.45X10 -8 m 2 /N (l.OxlO -4 (psi) -1 ) (ref. 12) 

The spinning torque for the naphthenic mineral oil was similarly examined by both 
theories. The parameters which gave the best fit are 

a = 4.93X10 -8 m 2 /N (3.4X10 -4 (psi) -1 ) 
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0 50 100 150 

Pressure, ksi 


Figure 5. - Pressure-viscosity model for di -2 -ethyl hexyl sebacate obtained 
from spinning torque data taken in spinning torque apparatus. Test con- 
ditions: temperature, ambient; spinning speed, 950 rpm ; ball -groove 
conformity 51 percent; ambient viscosity, p 0 , 16xl0" 3 N-secfm 2 {2.3xl0~ 6 
lb-sec/in. ); stress at which pressure -viscosity exponent changes, Si, 
275xl0 6 N/m 2 (40 000 psi). 


0 = o 


p = 0.079 N-sec/m 2 (1. 13xl0 -5 lb-sec/in. 2 ) 


For the conventional EHD analysis: 

Sj = 124X10 6 N/m 2 (18 000 psi) 


For the microasperity EHD analysis: 

= 131X10 6 N/m 2 (19 000 psi) 


11 




Maximum Hertz stress, N/m 2 


I I I I 1 ! I 1 L J 

60 80 100 120 140 60 80 100 120 140 

Maximum Hertz stress, ksi 

(a) Conventional elastohyd nodynamic analysis. (b) Microasperity el a sto hydrodynamic analysis. 

Figure 6. - Comparison of experimental and theoretical spinning torque for di-2-ethylhexyl sebacate. Pressure -viscosity coefficients: a, 2.03x10 $ m 2 /l\I 
(1.4x10^ (psiP 1 ); |3, 1.45xl0 -10 m 2 /N (1. lxlO' 6 (psiP 1 ). Ambient viscosity, p^, 16x10^ N-sec/m? (2.3xl0~ 6 Ib-sec/in . 2 ); stress at which pressure- 
viscosity exponent changes, 275xl0 6 N/m 2 (40 000 psil. 


The preceding values yield the lubricant model shown in figure 7. The spinning torque 
for the two analyses is shown in figure 8. The value of a used here is considerably 
higher than the previously published values 

a = 0.94X10 -8 m 2 /N (0.65X10 -4 (psi)' 1 ) (ref. 11) 

and 

a- = 2.15X10 -8 m 2 /N (1.48X10" 4 (psi) -1 ) (ref. 13) 
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0 50 100 150 

Pressure, ksi 

Figure 7. - Pressure-viscosity models for super-refined napthenic min- 
eral oil obtained from spinning torque data taken in spinning torque 
apparatus. Test conditions: temperature, ambient; spinning speed, 

950 rpm; ball -groove conformity, 51 percent. 


Although the same value of a was used for the conventional and microasperity analyses, 
it was necessary to use different values of Sj for the two analyses to obtain the best 
agreement. 

The spinning torque for the polyphenyl ether (5P4E) was also examined by the use of 
both analyses. The best agreement was obtained by using the following lubricant param- 
eters: 


a = 5.07X10" 8 m 2 /N (3.5X10 -4 (psi)" 1 ) 

H o = 0.8 N-sec/m 2 (1.16X10 -4 lb-sec/in. 2 ) 
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Maximum Hertz stress, ksl 

(a) Conventional elastohydrodynamic analysis; stress (b) Microasperity elastohydrodynamic analysis; stress 

at which pressure -viscosity exponent changes, S j, at which pressure -viscosity exponent changes, Sj, 

124xl0 6 N/m 2 (18 000 psi). 131X10 6 N/m 2 (19 000 psi). 

Figure 8. - Comparison of experimental and theoretical spinning torques for a super-refined napthenic mineral oil. Pressure -viscosity 
coefficients: a, 4.93x10"^ m^/N (3.4x10^ (psi)’l); p = 0. Ambient viscosity, 0.079 N-sec/m^ (l. 13x10"^ lb -sec/in. h. 

For the conventional EHD analysis: 

j3 = -1.16X10" 9 m 2 /N (-8.0X10 -6 (psi) -1 ) 

Sj = 117X10 6 N/m 2 (20 000 psi) 

For the microasperity EHD analysis: 

j3 = -1.4X10 -12 m 2 /N (-1.0X10 -8 (psi) -1 ) 

= 138X10 6 N/m 2 (20 000 psi) 

14 
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In this case, reasonable agreement can only be obtained by using different values of f 3 
and S 1 for the two different lubrication analyses. Also the value of 0 is negative in 
both cases, indicating an apparent reduction in viscosity at high pressures. The restat- 
ing lubricant pressure -viscosity models are portrayed in figure 9. The spinning torques 



0 50 100 150 

Pressure, ksi 


Figure 9. - Pressure -viscosity models used for polypheny I ether (5P4E) 
obtained from spinning torque data taken in spinning torque apparatus. 
Test conditions: temperature, ambient; spinning speed, 930 rpm ; ball- 
groove conformity, 51 percent. 
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60 80 100 120 60 80 100 120 
Maximum Hertz stress, ksi 

(a) Conventional elastohydrodynamic analysis; (b) Microasperity elastohydrodynamic analysis; 
pressure -viscosity coefficient (3, -1.16x10-9 pressure -viscosity coefficient, [3, -1.4xl0~12 
m2/N (-8. 0x10-6 (psi)-l); stress at which m 2 /N (-1. Ox 10 -8 (psi)“l); stress at which 

pressure-viscosity exponent changes, Sj, pressure-viscosity exponent changes, S|, 

117X10 6 N/m 2 (17 000 psi). 138xl0 6 N/m 2 (20 000 psi). 

Figure 10. - Comparison of theoretical and experimental spinning torques for polyphenyl ether (5P4E). Pressure- 
viscosity coefficient o, 5.07x10"^ m 2 /N (3.5x10"^ (psiP^); ambient viscosity. Mo- 0-8 N-sec/m 2 (1. 16x10 
Ib-sec/in. 2 ). 


against maximum Hertz stress for the two lubrication analyses are shown in figure 10. 

The value of a used in this analysis again differs from the previously published 
values 

a = 2xl0 -8 m 2 /N (1.38X10 -4 (psi)" 1 ) (ref. 11) 

and 

a = 4.64X10 -8 m 2 /N (3.2xl0~ 4 (psi) -1 ) (ref. 14) 

The explanation probably lies in the extremely large shear rates present in the spinning 
ball system. 
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For comparison purposes the spinning torque theory for the microasperity EHD 
analysis was used with the lubricant model for the conventional EHD analysis to calcu- 
late the spinning torque for both the super -refined naphthenic mineral oil and the poly- 
phenyl ether. The results of these comparisons are presented in figure 11. From fig- 
ure 11(a) it can be seen that the calculated spinning torque for the super -refined naph- 



60 80 100 120 60 80 100 120 
Maximum Hertz stress, ksi 

(a) Super-refined napthemc mineral oil. Pressure- (b) Polyphenyl ether (5P4E). Pressure -viscosity 
viscosity coefficients: a, 4.93xl0~8 m 2 /l\| (3.4xl(H coefficients: a, 5. 07x10"^ m 2 /N (3. 5x10'^ fpsirt; 
(psir^i; (5, 0. Ambient viscosity, pg, 0. 074 N- sec/ p, -1. 16x10*9 m 2 /N (-8. 0xl0'6 ipsirt. Ambient 

m 2 (1. nxlO'-Mb-sec/in. 2 i; stress at which viscosity, pg, 0.8 N-sec/m 2 (1. 16 xKH Ib-sec/in. ?); 

pressure -viscosity exponent changes, Si, 124x10^ stress at which pressure -viscosity exponent 
N/m 2 (18 000 psi). changes, S\, 117x10-6 N/m 2 (17 000 psi). 


Figure 11. - Comparison of experimental and theoretical spinning torques using a microasperity analysis and same lubricant 
model as for conventional elastohydrodynamic lubrication. 


thenic mineral oil is somewhat less than the experimental spinning torque. In fig- 
ure 11(b) the calculated spinning torque for the polyphenyl ether is much less than the 
experimental spinning torque. The basic reason for the lower calculated spinning 
torques when using this method is the difference in EHD film thickness between the 
microasperity and conventional EHD analyses. The film thickness was calculated for all 
three lubricants by the method of equation (3) for the conventional EHD analysis and by 
equation (10) for the microasperity EHD analysis. The results are plotted in figure 12, 
from which it is readily seen that the conventional EHD film thickness is roughly 
50 times the microasperity EHD film thickness. 
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Film thickness, 



CONCLUDING REMARKS 


The results presented herein have demonstrated the possibility of either conventional 
elastohydrodynamic or microasperity elastohydrodynamic lubrication for the various 
lubricants tested. The values of pres sure -viscosity coefficient which produce the best 
fit between the experimental data and the analysis differ somewhat from the published 
values. The published values for the initial pres sure -viscosity exponent are by no 
means consistent. The differences between them are as great or greater than the differ- 
ence between the present values and published values. 

The apparent change inthe pressure -viscosity exponent (from a to /3) must occur 
regardless of which EHD analysis is used. Without a reduction in this exponent the 
torque values become impossibly high. This reduction has been rarely observed in the 
past, probably because most of the previous traction data were taken at much lower 
Hertz stresses and lower shear rates. A somewhat similar change in the exponent is, 
however, given in reference 15, although the transition pressure in reference 15 is much 
higher . 

The polyphenyl ether exhibits a peculiar characteristic in that a negative secondary 
slope /3 is required for the pressure-viscosity curve to give reasonable correlation be- 
tween the theoretical and experimental results. In other rolling-element bearing studies 
(refs. 16 to 18), polyphenyl ether has proved to be a poor lubricant at high Hertz 
stresses. This may be due to the apparent reduction in viscosity observed in the cur- 
rent work. 

The theoretical work reported herein has been on the basis of either a conventional 
elastohydrodynamic analysis or a microasperity elastohydrodynamic analysis with no 
coupling between the two. In a real system there must be an interaction of the conven- 
tional and microasperity elastohydrodynamic lubrication effects. The probable result of 
this would be to increase the film thickness. However, in a real system the side leakage 
effects would tend to reduce the film thickness. The effect of neglecting coupling be- 
tween the two systems and side leakage are thus opposite. 


SUMMARY OF RESULTS 

The microasperity analysis for elastohydrodynamic lubrication of a ball spinning on a 
flat surface was extended to a ball spinning in a nonconforming groove. The analysis was 
compared to the conventional elastohydrodynamic analysis of a ball spinning in a noncon- 
forming groove. Using both analyses, a rheological model was constructed for a di-2- 
ethylhexyl sebacate, a super -refined naphthenic mineral oil, and a polyphenyl ether 
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(5P4E) based upon spinning torque data obtained in another investigation. The following 
results were obtained: 

1. The value of the lubricant pressure -viscosity coefficient that made the data fit the 
analysis of the fluids differs somewhat from published values. For all three lubricants, 
an exponential composite model, rather than a single exponential model, best described 
the lubricant rheology . 

2. A reasonably good agreement existed between the experimental and analytical 
values of torque when either the conventional or the microasperity elastohydrodynamic 
analysis for a ball spinning in a nonconforming groove was used. 

3. The polyphenyl ether lubricant exhibited a negative secondary slope in the com- 
posite exponential model, which would imply a decrease in viscosity with pressure at 
pressures beyond 117x10® N/m^ (17 000 psi). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 26, 1972, 

132-15. 
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APPENDIX - SYMBOLS 


2 2 

A total asperity area, m (in. ) 

a major semiaxis of contact ellipse, m (in.) 

b minor semiaxis of contact ellipse, m (in.) 

b' semiwidth of contact ellipse at y, m (in.) 

9 

E modulus of elasticity, N/m (psi) 

E’ materials properties factor, N/m (psi) 

F friction force 

h film thickness, m (in.) 

h Q minimum distance between ball and groove, m (in.) 

K constant defining outer boundary of integration 

M torque due to viscous drag, N-m (lb-in.) 

M g total spinning torque in Hertzian ellipse, N-m (lb -in.) 

Mj spinning torque in Hertzian ellipse outside inscribed circle, N-m (lb-in.) 

Mg spinning torque in Hertzian ellipse inside inscribed circle, N-m (lb-in.) 

P normal load, N (lb) 

2 

average pressure over asperity, N/m (psi) 
p pressure, N/m (psi) 

R radius of equivalent cylinder, m (in.) 

Rg radius of ball, m (in.) 

Rq radius of groove, m (in.) 

r ,<p polar coordinates 

r Q value of r at outer boundary of Hertzian ellipse, m (in.) 

o 

S contact stress, N/m (psi) 

Sj stress at which pressure -viscosity exponent changes, N/m 2 (psi) 

v entraining velocity, m/sec (in. /sec) 

W load per unit width, N/m (lb /in.) 

x,y,z Cartesian coordinates, m (in.) 
a, (3 pressure -viscosity exponents 
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a Poisson's ratio 

p 

T shear stress, N/m (psi) 

2 2 

n absolute viscosity, N-sec/m (lb-sec/in. ) 

p equivalent viscosity, N-sec/m 2 (lb -sec /in. 2 ) 

p p 

p 0 ambient viscosity , N-sec/m (lb-sec/in. ) 

P b viscosity at y = b, N-sec/m (lb-sec/in. 2 ) 

p radius of curvature of asperity, m (in.) 

co angular velocity, rad/sec 
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